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INTRODUCTION 

A  set  of  subroutines  has  been  developed  at  the  AEC 
Computing  Facility  at  Nev;  York  University  for  evaluating 
simple  functions  and  for  performing  certain  s-^andard  pro- 
cedures.  Detailed  descriptions  of  the  routines  are  avail- 
able, and  a  library  tape  containing  the  routines  Is  kept 
at  the  Unlvac.  ^i!he   routines  are  written  to  be  compatible 
with  the  NYU  Compiler  System  (NYO-6478).   The  num.ber  of 
routines  is  contantly  being  Increased  and  the  routines 
themselves  may  be  m.odlfled  from  time  to  tim.e  on  the  basis 
of  experience  in  applying  thorn  to  various  problems.   The 
llfst  presented  here  includes  those  routined  currently  in 
use,  together  with  a  brief  description  of  each.   A  revised 
index  of  current  routines  will  be  issued  at  intervals,  but 
before  using  any  of  the  routines  the  programmer  should  check 
the  library  subroutine  book  near  the  computer  for  the  latest 
corrections  and  modifications. 
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F.D.C.  Indicates  floating  decimal  composed  arithmetic. 
F.D.D.     "        "        "    :.D  composed    " 
P.P.       "     fixed  point  arithmetic. 


ARITHMETIC  ROUTINES 

Interpretive  Floating  £ecim.al  Routine  P-550106 

(AFOOI6  00000;3^) 

A  routine  to  interpret  successive  lines  of 
three-address  Injtruotions  ao  composed  float- 
ing decimal  operations. 

Composed  Floating  Decimal  P-55IIO8 

(CDOO60  00000/) 

A  one-block  routine  for  performing  addition, 

multiplication,  division  and  composition  of 

numbers  expressed  in  composed  floating  decimal 

form. 

Floating  Point  Arithmetic  P-5505II 

(FP0120  MAR55?^) 

A  routine  similar  to  th^  commonly  used  float- 
ing decimal  routine,  FR,  save  that  it  is  speeded 
up  by  eliminating  intermediate  composition  of 
numbers,  and  is  shorted  to  the  extent  that  fewer 
R  instructions  need  bo  coded. 

Floating  Decimal  Routine  P-5504o4 

(FR0120  00000/) 

This  is  the  floating  point  routine  which  has  been 

in  use  the  longest  at  NY^J.   The  designations 

F.D.C.  (Floating  Deci:r.ul  Composed)  and  F.D.D. 

(Floating  Decimal  Decomposed)  in  pax-'enthes^s 

after  other  subroutines  imply  that  FR  must  be 

compiled  into  lines  830  snd  following. 


-h. 


Complex  Nu-T-ocr  Arithmetic  (F.D.C.)  P- 551 10? 

(CEOO60  00000/; 

Addition,  7'ubtrc. oticn,  multiplication,  and 

division  of  ccniple?:  nunberD. 

Complex  NiTaber  Aritl'imstic  (?.?.;  P-551106 

(CXOOoO  00C00;L0 

Addition,  multiplication,  and  division  of 

c  omp 1 9 X  numb e  r  s . 

BLOCK  0r~PUT  R0UTIN3S 

Block  Output  P-550109 

(BLOOIO  BLOCK;^; 

On  each  entrance  to  this  routine  the  quantity 

in  rA  is  stored  at  ':h3  next  available  location 

in  a  specified  output  block.   When  the  output 

block  is  filled  it  is  v;ritten  on  a  designated 

tape. 

Block  Output  (V-register^  P-55OIO8 

(BVOOlO  BLOCK/) 

A  roucine  similar  to  BL  except  that  the  block 

is  filled  by  pairs  of  quantities  from  rV. 

DEBUGGING  ROUTIKES 

Generator  for  Debugging  Tape  P- 560104 

(GROI80  00000/) 

A  routine  to  generate  a  requested  debugging 

routine  into  available  blocks  of  storage  in 

a  subject  program. 
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Dsbuggli-r;  F.cutir.e  P-55051C 

(DEOC6O  OOOCO,t') 

A  rcuvj.jf.e  to  porxcrm  automatic  storage  dumpc 

(17  blc^kr.  cr.to  ""apj:;  or  to  s-ve  automatic 

i'-eglGter  ard  st-or-yg^  type-outs. 

Storage  Dur::p  T'cut-.n3S 

(DS0008  ccoo  ;■ 

(DPC026  00000^') 

DS  dUTips  17  blocko  o::'  storso-s  onto   servo  A. 

DP  preserves  the  contents  of  registers,  dumps 

a  selected  region  of  storage  onto  a  selected 

servo,  and  restores  the  registers. 

Q,T,U  Tracer  -  SC  typevjrlter  output  P-56OIO5 

(TROO6O  00000/) 

A  routine  to  rr.onlto:.^  a  subject  program, 

recording  che  transfer  Instructions  as  they 

are  exeouued. 

Q,T,U  Tracer  -  tape  output  P-560IO6 

(TEO229  00000/) 

A  routine  similar  to  the  above  TR. 


Trac er  Counter 

(TC0120  ooo.yj/) 

A  routine  to  monitor  a  subject  program, 

recording  the  changes  effected  on  certain 

given  line  locations. 

rA  Monitor  -  SC  typevJi-iter  output 
(KN0102  OCCOO/j 

A  routine  to  monitor  a  subject  program, 
recording  the  changes  ofioctoo  on  the 
A  register. 


P-560107 


P-560108 
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rA  Monitor  -  tape  output  P-56OIO9 

(MEOII9  00000^) 

A  routine  similar  to  the  above  MN. 


Follower 

(FL0120  00000^) 

A  compilable  version  of  the  NYU  service  routine 

called  Follower. 


P-5^0729-2 


DIFFERENTIAL  EQUATION  ROUTINF^ 

Differential  Equation  Solution  (F.P.)  P-5*^08lO 

(DYOO62  GILLO^) 

Integration  of  a  system  of  ordinary  differen- 
tial equations  using  the  Runge-Kutta-Glll 
method . 

Differential  Equation  Solution  P-550512 

(DFOO69  FLOATS) 

Integration  of  a  system  of  ordinary  differen- 
tial equations  using  the  Runge-Kutta-Glll 
method.  The  arithmetic  is  floating  decimal 
requiring  use  of  routine  F.P. 

EDIT  ROUTINES 

Sign  Edit  P-540830 

(EAOOO7  +SIGN^) 

A  routine  to  prefix  a  positive  number  with  a 

plus  sign  in  place  of  the  Univac  zero. 

Edit  Floating  Decimal  Composed  Number  P-55010^ 

(EBOO34  EDITA^) 

A  routine  to  edit  a  composed  floating  point 

number  into  two  words--a  signed  mantissa  with 

decimal  point  and  a  signed  base  10  exponent. 
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Edit  Floating  Decimal  Decomposed  Nvjnber 
(EPO035  EDITA/) 

A  routine  to  edit  a  decomposed  floating  point 
n^ambor  into  a  signed  mantissa  with  decimal 
point  and  a  signed  base  10  exponent. 


P-550110 


EXPONENTIALS  AND  ROOTS 

A^  (7.?. ) 

(AXOO32  00000/) 

A  routine  to  evaluate  A^  for 
0  <  A  <  1 
0  <  ::  <  1  . 


P-550414 


,x 


10 


X 


P-5504l6 


(LH0040  EXP 10/) 

A  routine  to  evalua-:e  10"^  for  |xj  <  50  with 


^k 


error  j-.  COOO8  10  (k  =  fractional  part  of  |x|). 

e^  (F.D.D. ) 
(LPOO52  EXPOO/) 


e^  (?.P.)(x  non-positive) 
(LTOO55   .POO/) 

Square  Pool  (F.D.C. ) 
(NA002iL  R0Gr2/) 

Square  Root  (F.D.D) 
(NJOO23  R00T2/) 

Square  Root  (F.P. ) 
(NS0012  R0CT2/) 

Integral  Root;  A  ^^ 

(KX003^  ApA:.:.T/) 

|A|  <  1;  p  >  2 

(Newton-Raphscn  iterative  procedure.) 
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FILL  ROUTINES 


Line  Fill  Routlna  P-54o6o4 

(F()000()  AFILL;J^) 

A  set  of  routines  which  can  be  used  In  com- 
piling to  fill  from  one  to  nine  lines  with 
zeros. 

IGNORE  SQUEEZE  ROUTINES 

Ignore  Squeezo  Block  Output  P-550329 

(GN0057  00000^) 

On  each  entrance  to  this  routine  the  ignores 

are  removed  from  the  word  in  rA  and  the 

result  is  stored  into  an  output  block. 

(Faster  than  routine  IG. ) 

Ignore  Squeeze  Block  Output  P-551105 

(IG003^  00000^) 

Same  as  routine  GN  in  purpose  but  shorter 

in  space. 

INTEGRATION  ROUTINES 

Integration:    Sim.pson's  Rule    (F.D.D.)  P-5^0607 

(QNOO5O  00000^)  r\> 

A  routine  to  evaluate  \    f(x)dx 

y  a 

given  a,  b,  the  interval  of  Integration 

and  a  routine  for  f(x). 

Integration:   Gnuss's  Method  (F.D.D.)  P-550506 

(GMOO36  00000^) 

Evaluation  of  an  integral  using  Gauss's 

method  (integrand  evaluated  at  unequal 

intervals) of  order  n  (4  <  n  <  lO). 
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Integration:  Gauss's  Method  (P.P.)  P-550505 

(GS0025  00000^) 

Fixed  point  version  of  the  above  routine, 

GM. 

G4  * • *  through  G-  are  the  routines  containing 
the  coefficients  for  GM  and  GS. 

INTERPOLATION  ROUtlt^ES 


Interpolation  (F.D.D.)  P-54o8ll 

(IPOO62  AITK.;i^) 

A  routine  to  perform  interpolation  of  order 

n  using  Aitkeh's  method  (resti*icted  to  equal 

intervals) . 

Lagharige  Intefpolation  (F.D.D.)  P-540624 

(ILOO69  00000^) 

A  routine  to  perform  interpolation  of  order 

n  using  Lagrange's  formula  (untecjual  intervals). 

LOGARITHM  ROUTINES 

Natural  Logarithm  (F-D.C.)  P-550830 

(LF0044  LNOOOX) 

Logarithm  to  Base  10  (p'^'d*)  P-550415 

(LN0022  LOGlOjif) 

A  routine  to  find  log,QX  with  | error |  <  .0007- 

Logarithm  to  Base  2  (F.D.C.)  P-540702 

(LJ0040  BASE2;lf) 

Logarithm  to  Base  2  (P.P.)  P-5^0712 

(LX0043  BASE  2X) 
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POLYNOMIAL  ROUTINES 

Polynomial  Evaluation  (F.P.) 
(PS0018  00000^)  ^ 

A  routine  to  evaluate 

n,  the  a  . ' s  and  x. 
J 

SORTING  ROUTINES 

Sorting  Routine 
(ST0032  SORTA/) 


J=0 


ajx^ 
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given 


A  routine  to  sort  Unlvac  words  (numbers 
and/or  letters)  into  a  sequence  of  increasing 
magnitude.   Up  to  968  words  may  be  sorted. 


P-5if0606 


P-550827 


TRIGONOMETRIC  ROUTINES 

Sine-Cosine  (F.D.C.) 

(TDOO52  COSSN;J^) 

Sine-Cosine  (F.D.D.) 
(TKOO54  COSSNiJ^) 

Sine-Cosine  (F.P.) 
(TV0041  00000;j^) 

Tan"  X  (F.P. ) 

(TY0057  OOCOO^) 

Continued  fraction  approximation,  error  <  10 


P-55OIO7 


P-55OIO5 


P-550413 


P-550902 


-10 


Tan~  X  (F.P. ) 

(TZOO25  00000^) 

Polynomial  approximation,  error  ^  10 

(Faster  but  less  accurate  than  TY) 


P-550901 
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